Perspective on Integrative Simulations of Bioenergetic Domains.

Bioenergetic processes in cells, such as photosynthesis or respiration, integrate many time and length scales, which makes the simulation of energy conversion with a mere single level of theory impossible. Just like the myriad of experimental techniques required to examine each level of organization, an array of overlapping computational techniques is necessary to model energy conversion. Here, a perspective is presented on recent efforts for modeling bioenergetic phenomena with a focus on molecular dynamics simulations and its variants as a primary method. An overview of the various classical, quantum mechanical, enhanced sampling, coarse-grained, Brownian dynamics, and Monte Carlo methods is presented. Example applications discussed include multiscale simulations of membrane-wide electron transport, rate kinetics of ATP turnover from electrochemical gradients, and finally, integrative modeling of the chromatophore, a photosynthetic pseudo-organelle.

The electrostatic landscape of MHC-peptide binding revealed using inception networks.

Predictive modeling of macromolecular recognition and protein-protein complementarity represents one of the cornerstones of biophysical sciences. However, such models are often hindered by the combinatorial complexity of interactions at the molecular interfaces. Exemplary of this problem is peptide presentation by the highly polymorphic major histocompatibility complex class I (MHC-I) molecule, a principal component of immune recognition. We developed human leukocyte antigen (HLA)-Inception, a deep biophysical convolutional neural network, which integrates molecular electrostatics to capture non-bonded interactions for predicting peptide binding motifs across 5,821 MHC-I alleles. These predictions of generated motifs correlate strongly with experimental peptide binding and presentation data. Beyond molecular interactions, the study demonstrates the application of predicted motifs in analyzing MHC-I allele associations with HIV disease progression and patient response to immune checkpoint inhibitors. A record of this paper's transparent peer review process is included in the supplemental information.

Structural mechanism of HP1α-dependent transcriptional repression and chromatin compaction.

Heterochromatin protein 1 (HP1) plays a central role in establishing and maintaining constitutive heterochromatin. However, the mechanisms underlying HP1-nucleosome interactions and their contributions to heterochromatin functions remain elusive. In this study, we employed a multidisciplinary approach to unravel the interactions between human HP1α and nucleosomes. We have elucidated the cryo-EM structure of an HP1α dimer bound to an H2A.Z nucleosome, revealing that the HP1α dimer interfaces with nucleosomes at two distinct sites. The primary binding site is located at the N-terminus of histone H3, specifically at the trimethylated K9 (K9me3) region, while a novel secondary binding site is situated near histone H2B, close to nucleosome superhelical location 4 (SHL4). Our biochemical data further demonstrates that HP1α binding influences the dynamics of DNA on the nucleosome. It promotes DNA unwrapping near the nucleosome entry and exit sites while concurrently restricting DNA accessibility in the vicinity of SHL4. This study offers a model that explains how HP1α functions in heterochromatin maintenance and gene silencing, particularly in the context of H3K9me-dependent mechanisms. Additionally, it sheds light on the H3K9me-independent role of HP1 in responding to DNA damage.

Structural insights into the modulation of coronavirus spike tilting and infectivity by hinge glycans.

Coronavirus spike glycoproteins presented on the virion surface mediate receptor binding, and membrane fusion during virus entry and constitute the primary target for vaccine and drug development. How the structure dynamics of the full-length spikes incorporated in viral lipid envelope correlates with the virus infectivity remains poorly understood. Here we present structures and distributions of native spike conformations on vitrified human coronavirus NL63 (HCoV-NL63) virions without chemical fixation by cryogenic electron tomography (cryoET) and subtomogram averaging, along with site-specific glycan composition and occupancy determined by mass spectrometry. The higher oligomannose glycan shield on HCoV-NL63 spikes than on SARS-CoV-2 spikes correlates with stronger immune evasion of HCoV-NL63. Incorporation of cryoET-derived native spike conformations into all-atom molecular dynamic simulations elucidate the conformational landscape of the glycosylated, full-length spike that reveals a role of hinge glycans in modulating spike bending. We show that glycosylation at N1242 at the upper portion of the stalk is responsible for the extensive orientational freedom of the spike crown. Subsequent infectivity assays implicated involvement of N1242-glyan in virus entry. Our results suggest a potential therapeutic target site for HCoV-NL63.

Adaptive Ensemble Refinement of Protein Structures in High Resolution Electron Microscopy Density Maps with Radical Augmented Molecular Dynamics Flexible Fitting.

Recent advances in cryo-electron microscopy (cryo-EM) have enabled modeling macromolecular complexes that are essential components of the cellular machinery. The density maps derived from cryo-EM experiments are often integrated with manual, knowledge-driven or artificial intelligence-driven and physics-guided computational methods to build, fit, and refine molecular structures. Going beyond a single stationary-structure determination scheme, it is becoming more common to interpret the experimental data with an ensemble of models that contributes to an average observation. Hence, there is a need to decide on the quality of an ensemble of protein structures on-the-fly while refining them against the density maps. We introduce such an adaptive decision-making scheme during the molecular dynamics flexible fitting (MDFF) of biomolecules. Using RADICAL-Cybertools, the new RADICAL augmented MDFF implementation (R-MDFF) is examined in high-performance computing environments for refinement of two prototypical protein systems, adenylate kinase and carbon monoxide dehydrogenase. For these test cases, use of multiple replicas in flexible fitting with adaptive decision making in R-MDFF improves the overall correlation to the density by 40% relative to the refinements of the brute-force MDFF. The improvements are particularly significant at high, 2-3 Å map resolutions. More importantly, the ensemble model captures key features of biologically relevant molecular dynamics that are inaccessible to a single-model interpretation. Finally, the pipeline is applicable to systems of growing sizes, which is demonstrated using ensemble refinement of capsid proteins from the chimpanzee adenovirus. The overhead for decision making remains low and robust to computing environments. The software is publicly available on GitHub and includes a short user guide to install R-MDFF on different computing environments, from local Linux-based workstations to high-performance computing environments.

Functional Organization of Glycolytic Metabolon on Mitochondria.

Glucose, the primary cellular energy source, is metabolized through glycolysis initiated by the rate-limiting enzyme Hexokinase (HK). In energy-demanding tissues like the brain, HK1 is the dominant isoform, primarily localized on mitochondria, crucial for efficient glycolysis-oxidative phosphorylation coupling and optimal energy generation. This study unveils a unique mechanism regulating HK1 activity, glycolysis, and the dynamics of mitochondrial coupling, mediated by the metabolic sensor enzyme O-GlcNAc transferase (OGT). OGT catalyzes reversible O-GlcNAcylation, a post-translational modification, influenced by glucose flux. Elevated OGT activity induces dynamic O-GlcNAcylation of HK1's regulatory domain, subsequently promoting the assembly of the glycolytic metabolon on the outer mitochondrial membrane. This modification enhances HK1's mitochondrial association, orchestrating glycolytic and mitochondrial ATP production. Mutations in HK1's O-GlcNAcylation site reduce ATP generation, affecting synaptic functions in neurons. The study uncovers a novel pathway that bridges neuronal metabolism and mitochondrial function via OGT and the formation of the glycolytic metabolon, offering new prospects for tackling metabolic and neurological disorders.

Multidimensional quantitative phenotypic and molecular analysis reveals neomorphic behaviors of p53 missense mutants.

Mutations in the TP53 tumor suppressor gene occur in >80% of the triple-negative or basal-like breast cancer. To test whether neomorphic functions of specific TP53 missense mutations contribute to phenotypic heterogeneity, we characterized phenotypes of non-transformed MCF10A-derived cell lines expressing the ten most common missense mutant p53 proteins and observed a wide spectrum of phenotypic changes in cell survival, resistance to apoptosis and anoikis, cell migration, invasion and 3D mammosphere architecture. The p53 mutants R248W, R273C, R248Q, and Y220C are the most aggressive while G245S and Y234C are the least, which correlates with survival rates of basal-like breast cancer patients. Interestingly, a crucial amino acid difference at one position-R273C vs. R273H-has drastic changes on cellular phenotype. RNA-Seq and ChIP-Seq analyses show distinct DNA binding properties of different p53 mutants, yielding heterogeneous transcriptomics profiles, and MD simulation provided structural basis of differential DNA binding of different p53 mutants. Integrative statistical and machine-learning-based pathway analysis on gene expression profiles with phenotype vectors across the mutant cell lines identifies quantitative association of multiple pathways including the Hippo/YAP/TAZ pathway with phenotypic aggressiveness. Further, comparative analyses of large transcriptomics datasets on breast cancer cell lines and tumors suggest that dysregulation of the Hippo/YAP/TAZ pathway plays a key role in driving the cellular phenotypes towards basal-like in the presence of more aggressive p53 mutants. Overall, our study describes distinct gain-of-function impacts on protein functions, transcriptional profiles, and cellular behaviors of different p53 missense mutants, which contribute to clinical phenotypic heterogeneity of triple-negative breast tumors.

High-resolution structures with bound Mn2+ and Cd2+ map the metal import pathway in an Nramp transporter.

Transporters of the Nramp (Natural resistance-associated macrophage protein) family import divalent transition metal ions into cells of most organisms. By supporting metal homeostasis, Nramps prevent diseases and disorders related to metal insufficiency or overload. Previous studies revealed that Nramps take on a LeuT fold and identified the metal-binding site. We present high-resolution structures of Deinococcus radiodurans (Dra)Nramp in three stable conformations of the transport cycle revealing that global conformational changes are supported by distinct coordination geometries of its physiological substrate, Mn2+, across conformations, and by conserved networks of polar residues lining the inner and outer gates. In addition, a high-resolution Cd2+-bound structure highlights differences in how Cd2+ and Mn2+ are coordinated by DraNramp. Complementary metal binding studies using isothermal titration calorimetry with a series of mutated DraNramp proteins indicate that the thermodynamic landscape for binding and transporting physiological metals like Mn2+ is different and more robust to perturbation than for transporting the toxic Cd2+ metal. Overall, the affinity measurements and high-resolution structural information on metal substrate binding provide a foundation for understanding the substrate selectivity of essential metal ion transporters like Nramps.

Integrated analyses reveal a hinge glycan regulates coronavirus spike tilting and virus infectivity.

Coronavirus spike glycoproteins presented on the virion surface mediate receptor binding, and membrane fusion during virus entry and constitute the primary target for vaccine and drug development. How the structure dynamics of the full-length spikes incorporated in viral lipid envelope correlates with the virus infectivity remains poorly understood. Here we present structures and distributions of native spike conformations on vitrified human coronavirus NL63 (HCoV-NL63) virions without chemical fixation by cryogenic electron tomography (cryoET) and subtomogram averaging, along with site-specific glycan composition and occupancy determined by mass spectroscopy. The higher oligomannose glycan shield on HCoV-NL63 spikes than on SARS-CoV-2 spikes correlates with stronger immune evasion of HCoV-NL63. Incorporation of cryoET-derived native spike conformations into all-atom molecular dynamic simulations elucidate the conformational landscape of the glycosylated, full-length spike that reveals a novel role of stalk glycans in modulating spike bending. We show that glycosylation at N1242 at the upper portion of the stalk is responsible for the extensive orientational freedom of the spike crown. Subsequent infectivity assays support the hypothesis that this glycan-dependent motion impacts virus entry. Our results suggest a potential therapeutic target site for HCoV-NL63.

Revealing a Hidden Intermediate of Rotatory Catalysis with X-ray Crystallography and Molecular Simulations.

The mechanism of rotatory catalysis in ATP-hydrolyzing molecular motors remains an unresolved puzzle in biological energy transfer. Notwithstanding the wealth of available biochemical and structural information inferred from years of experiments, knowledge on how the coupling between the chemical and mechanical steps within motors enforces directional rotatory movements remains fragmentary. Even more contentious is to pinpoint the rate-limiting step of a multistep rotation process. Here, using vacuolar or V1-type hexameric ATPase as an exemplary rotational motor, we present a model of the complete 4-step conformational cycle involved in rotatory catalysis. First, using X-ray crystallography, a new intermediate or "dwell" is identified, which enables the release of an inorganic phosphate (or Pi) after ATP hydrolysis. Using molecular dynamics simulations, this new dwell is placed in a sequence with three other crystal structures to derive a putative cyclic rotation path. Free-energy simulations are employed to estimate the rate of the hexameric protein transformations and delineate allosteric effects that allow new reactant ATP entry only after hydrolysis product exit. An analysis of transfer entropy brings to light how the side-chain-level interactions transcend into larger-scale reorganizations, highlighting the role of the ubiquitous arginine-finger residues in coupling chemical and mechanical information. An inspection of all known rates encompassing the 4-step rotation mechanism implicates the overcoming of the ADP interactions with V1-ATPase to be the rate-limiting step of motor action.

Characterization and computational simulation of human Syx, a RhoGEF implicated in glioblastoma.

Structural discovery of guanine nucleotide exchange factor (GEF) protein complexes is likely to become increasingly relevant with the development of new therapeutics targeting small GTPases and development of new classes of small molecules that inhibit protein-protein interactions. Syx (also known as PLEKHG5 in humans) is a RhoA GEF implicated in the pathology of glioblastoma (GBM). Here we investigated protein expression and purification of ten different human Syx constructs and performed biophysical characterizations and computational studies that provide insights into why expression of this protein was previously intractable. We show that human Syx can be expressed and isolated and Syx is folded as observed by circular dichroism (CD) spectroscopy and actively binds to RhoA as determined by co-elution during size exclusion chromatography (SEC). This characterization may provide critical insights into the expression and purification of other recalcitrant members of the large class of oncogenic-Diffuse B-cell lymphoma (Dbl) homology GEF proteins. In addition, we performed detailed homology modeling and molecular dynamics simulations on the surface of a physiologically realistic membrane. These simulations reveal novel insights into GEF activity and allosteric modulation by the plekstrin homology (PH) domain. These newly revealed interactions between the GEF PH domain and the membrane embedded region of RhoA support previously unexplained experimental findings regarding the allosteric effects of the PH domain from numerous activity studies of Dbl homology GEF proteins. This work establishes new hypotheses for structural interactivity and allosteric signal modulation in Dbl homology RhoGEFs.

Energy landscape of the SARS-CoV-2 reveals extensive conformational heterogeneity.

Cryo-electron microscopy (cryo-EM) has produced a number of structural models of the SARS-CoV-2 spike, already prompting biomedical outcomes. However, these reported models and their associated electrostatic potential maps represent an unknown admixture of conformations stemming from the underlying energy landscape of the spike protein. As with any protein, some of the spike's conformational motions are expected to be biophysically relevant, but cannot be interpreted only by static models. Using experimental cryo-EM images, we present the energy landscape of the glycosylated spike protein, and identify the diversity of low-energy conformations in the vicinity of its open (so called 1RBD-up) state. The resulting atomic refinement reveal global and local molecular rearrangements that cannot be inferred from an average 1RBD-up cryo-EM model. Here we report varied degrees of "openness" in global conformations of the 1RBD-up state, not revealed in the single-model interpretations of the density maps, together with conformations that overlap with the reported models. We discover how the glycan shield contributes to the stability of these low-energy conformations. Five out of six binding sites we analyzed, including those for engaging ACE2, therapeutic mini-proteins, linoleic acid, two different kinds of antibodies, switch conformations between their known apo- and holo-conformations, even when the global spike conformation is 1RBD-up. This apo-to-holo switching is reminiscent of a conformational preequilibrium. We found only one binding site, namely that of AB-C135 remains in apo state within all the sampled free energy-minimizing models, suggesting an induced fit mechanism for the docking of this antibody to the spike.

The ugly, bad, and good stories of large-scale biomolecular simulations.

Molecular modeling of large biomolecular assemblies exemplifies a disruptive area holding both promises and contentions. Propelled by peta and exascale computing, several simulation methodologies have now matured into user-friendly tools that are successfully employed for modeling viruses, membranous nano-constructs, and key pieces of the genetic machinery. We present three unifying biophysical themes that emanate from some of the most recent multi-million atom simulation endeavors. Despite connecting molecular changes with phenotypic outcomes, the quality measures of these simulations remain questionable. We discuss the existing and upcoming strategies for constructing representative ensembles of large systems, how new computing technologies will boost this area, and make a point that integrative modeling guided by experimental data is the future of biomolecular computations.

Exploring cryo-electron microscopy with molecular dynamics.

Single particle analysis cryo-electron microscopy (EM) and molecular dynamics (MD) have been complimentary methods since cryo-EM was first applied to the field of structural biology. The relationship started by biasing structural models to fit low-resolution cryo-EM maps of large macromolecular complexes not amenable to crystallization. The connection between cryo-EM and MD evolved as cryo-EM maps improved in resolution, allowing advanced sampling algorithms to simultaneously refine backbone and sidechains. Moving beyond a single static snapshot, modern inferencing approaches integrate cryo-EM and MD to generate structural ensembles from cryo-EM map data or directly from the particle images themselves. We summarize the recent history of MD innovations in the area of cryo-EM modeling. The merits for the myriad of MD based cryo-EM modeling methods are discussed, as well as, the discoveries that were made possible by the integration of molecular modeling with cryo-EM. Lastly, current challenges and potential opportunities are reviewed.

Anionic Lipids Confine Cytochrome c2 to the Surface of Bioenergetic Membranes without Compromising Its Interaction with Redox Partners.

Cytochrome c2 (cyt. c2) is a major element in electron transfer between redox proteins in bioenergetic membranes. While the interaction between cyt. c2 and anionic lipids abundant in bioenergetic membranes has been reported, their effect on the shuttling activity of cyt. c2 remains elusive. Here, the effect of anionic lipids on the interaction and binding of cyt. c2 to the cytochrome bc1 complex (bc1) is investigated using a combination of molecular dynamics (MD) and Brownian dynamics (BD) simulations. MD is used to generate thermally accessible conformations of cyt. c2 and membrane-embedded bc1, which were subsequently used in multireplica BD simulations of diffusion of cyt. c2 from solution to bc1, in the presence of various lipids. We show that, counterintuitively, anionic lipids facilitate association of cyt. c2 with bc1 by localizing its diffusion to the membrane surface. The observed lipid-mediated bc1 association is further enhanced by the oxidized state of cyt. c2, in line with its physiological function. This lipid-mediated enhancement is salinity-dependent, and anionic lipids can disrupt cyt. c2-bc1 interaction at nonphysiological salt levels. Our data highlight the importance of the redox state of cyt. c2, the lipid composition of the chromatophore membrane, and the salinity of the chromatophore in regulating the efficiency of the electron shuttling process mediated by cyt. c2. The conclusions can be extrapolated to mitochondrial systems and processes, or any bioenergetic membrane, given the structural similarity between cyt. c2 and bc1 and their mitochondrial counterparts.

ChAdOx1 interacts with CAR and PF4 with implications for thrombosis with thrombocytopenia syndrome.

Vaccines derived from chimpanzee adenovirus Y25 (ChAdOx1), human adenovirus type 26 (HAdV-D26), and human adenovirus type 5 (HAdV-C5) are critical in combatting the severe acute respiratory coronavirus 2 (SARS-CoV-2) pandemic. As part of the largest vaccination campaign in history, ultrarare side effects not seen in phase 3 trials, including thrombosis with thrombocytopenia syndrome (TTS), a rare condition resembling heparin-induced thrombocytopenia (HIT), have been observed. This study demonstrates that all three adenoviruses deployed as vaccination vectors versus SARS-CoV-2 bind to platelet factor 4 (PF4), a protein implicated in the pathogenesis of HIT. We have determined the structure of the ChAdOx1 viral vector and used it in state-of-the-art computational simulations to demonstrate an electrostatic interaction mechanism with PF4, which was confirmed experimentally by surface plasmon resonance. These data confirm that PF4 is capable of forming stable complexes with clinically relevant adenoviruses, an important step in unraveling the mechanisms underlying TTS.

Poor Person's pH Simulation of Membrane Proteins.

pH conditions are central to the functioning of all biomolecules. However, implications of pH changes are nontrivial on a molecular scale. Though a rigorous microscopic definition of pH exists, its implementation in classical molecular dynamics (MD) simulations is cumbersome, and more so in large integral membrane systems. In this chapter, an integrative pipeline is described that combines Multi-Conformation Continuum Electrostatics (MCCE) computations with MD simulations to capture the effect of transient protonation states on the coupled conformational changes in transmembrane proteins. The core methodologies are explained, and all the software required to set up this pipeline are outlined with their key parameters. All associated analyses of structure and function are provided using two case studies, namely those of bioenergetic complexes: NADH dehydrogenase (complex I) and Vo domain of V-type ATPase. The hybrid MCCE-MD pipeline has allowed the discovery of hydrogen bond networks, ligand binding pathways, and disease-causing mutations.

Energy Landscape of the SARS-CoV-2 Reveals Extensive Conformational Heterogeneity.

Cryo-electron microscopy (cryo-EM) has produced a number of structural models of the SARS-CoV-2 spike, already prompting biomedical outcomes. However, these reported models and their associated electrostatic potential maps represent an unknown admixture of conformations stemming from the underlying energy landscape of the spike protein. As for any protein, some of the spike's conformational motions are expected to be biophysically relevant, but cannot be interpreted only by static models. Using experimental cryo-EM images, we present the energy landscape of the spike protein conformations, and identify molecular rearrangements along the most-likely conformational path in the vicinity of the open (so called 1RBD-up) state. The resulting global and local atomic refinements reveal larger movements than those expected by comparing the reported 1RBD-up and 1RBD-down cryo-EM models. Here we report greater degrees of "openness" in global conformations of the 1RBD-up state, not revealed in the single-model interpretations of the density maps, together with conformations that overlap with the reported models. We discover how the glycan shield contributes to the stability of these conformations along the minimum free-energy pathway. A local analysis of seven key binding pockets reveals that six out them, including those for engaging ACE2, therapeutic mini-proteins, linoleic acid, two different kinds of antibodies, and protein-glycan interaction sites, switch conformations between their known apo- and holo-conformations, even when the global spike conformation is 1RBD-up. This is reminiscent of a conformational pre-equilibrium. We found only one binding pocket, namely antibody AB-C135 to remain closed along the entire minimum free energy path, suggesting an induced fit mechanism for this enzyme.

Large-Scale Molecular Dynamics Simulations of Cellular Compartments.

Molecular dynamics or MD simulation is gradually maturing into a tool for constructing in vivo models of living cells in atomistic details. The feasibility of such models is bolstered by integrating the simulations with data from microscopic, tomographic and spectroscopic experiments on exascale supercomputers, facilitated by the use of deep learning technologies. Over time, MD simulation has evolved from tens of thousands of atoms to over 100 million atoms comprising an entire cell organelle, a photosynthetic chromatophore vesicle from a purple bacterium. In this chapter, we present a step-by-step outline for preparing, executing and analyzing such large-scale MD simulations of biological systems that are essential to life processes. All scripts are provided via GitHub.